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ABSTRACT 


Integral  equations  for  the  current  distribution  on  an  arbitrary 
wire  antenna  have  been  programmed  for  a  digital  computer.  The 
resulting  computer  program  will  calculate  the  current  distribution, 
the  input  impedance,  and  the  radiation  pattern  of  a  wire  antenna  of 
arbitrary  geometry.  Considered  as  a  transmitting  antenna,  the 
antenna  is  excited  by  electromotive  forces  applied  in  any  number  of 
gap i  along  the  wires.  The  program  developed  will  permit  these  gaps 
to  be  connected  by  an  arbitrary  non-radiating  network,  and  it  will 
take  into  account  resistive  and  reactive  loading  along  the  wireB.  The 
program  will  also  treat  the  antenna  as  a  receiving  antenna,  if  desired, 
and  then  uses  a  distant  dipole  as  the  source  of  the  induced  current 
distribution.  The  numerical  integral  equation  method  used  is  essen¬ 
tially  an  exact  method.  The  principal  errors  in  the  calculations 
arise  in  the  numerical  integration  of  the  integrals  in  the  integral 
equations,  and  these  errors  can  be  controlled.  It  is  usually  quite 
easy  to  obtain  accuracies  far  better  than  normal  measurement  ac¬ 
curacies,  and  at  a  much  lower  cost. 

Numerical  results  are  presented  for  a  simple  dipole,  a  Yagi 
antenna  array,  a  log-periodic  antenna  array,  an  airplane  tail-cap 
antenna,  a  center-fed  whip  antenna,  and  a  compact  VHF  antenna. 

Most  results  have  been  verified  by  measurements. 
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1. 


INTRODUCTION 


Most  theoretical  methods  available  today  for  calculating  the 
current  distribution,  the  input  impedance,  and  the  radiation  pattern 
of  wire  antennas  are  of  rather  limited  accuracy.  Most  methods  are 
not  concerned  with  an  evaluation  of  the  current  distribution  but  use, 
rather,  an  assumed  current  distribution  for  calculating  the  radiation 
pattern  and  the  input  impedance.  The  simple  electrical  dipole  is  the 
only  wire  antenna  for  which,  to  the  knowledge  of  the  authors,  attempts 
have  been  made  to  calculate  the  current  distribution  with  high  accuracy. 
Unfortunately,  however,  the  analytical  methods  used  to  calculate  the 
current  distribution  on  dipoles  are  not  generally  applicable  to  a  wire 
antenna  of  arbitrary  configuration.  The  reason  for  this  is  to  be  found 
in  the  limitations  of  the  so-called  analytical  methods.  For  many  years 
the  application  of  mathematical  analysis  has  been  the  only  method  that 
could  be  used  to  attempt  the  solution  of  electromagnetic  boundary- 
value  problems.  Today,  with  the  existence  of  fast  digital  computers, 
this  is  no  longer  true.  Numerous  problems  that  used  to  be  considered 
unsolvable,  or  for  which  only  approximate  solutions  could  be  obtained, 
can,  in  fact,  be  solved  today  in  great  generality  and  with  high  accuracy. 
The  general  wire  antenna  is  one  such  problem. 

Before  describing  the  method  presented  in  this  report  to  solve 
the  general  wire  antenna  problem,  we  shall  give  a  brief  account  of  the 
most  important  analytical  methods  that  have  been  used  to  calculate 
antenna  impedance.  The  best  known,  and  also  the  simplest,  analytical 
methods  for  calculating  the  input  impedance  of  a  wire  antenna  make 
use  of  an  assumed  distribution,  of  current  on  the  antenna.  By  the 
Poynting-vector  method,  one  calculates  the  radiated  power  by  inte¬ 
grating  the  far-zone  field  of  the  assumed  current  distribution.  The 
input  resistance  of  the  antenna  is  then  the  ratio  between  the  radiated 
power  and  the  square  of  the  input  current.  The  Poynting-vector  method 
was  used  by  Hertz*  as  early  as  1888  for  calculating  the  power  radiated 
by  a  Hertzian  dipole. 

Since  the  Poynting-vector  method  uses  only  the  far-zone  field 
for  the  impedance  calculation,  this  method  can  give  no  information  on 
the  antenna  reactance  which  is  associated  with  energy  stored  in  the  im¬ 
mediate  vicinity  of  the  antenna.  This  disadvantage  is  avoided  in  the  emf- 
method  which  uses  the  near-zone  field  for  the  impedance  calculation. 

The  emf-method  also  uses  an  assumed  current  distribution  on  the  antenna. 
This  current  distribution  radiates  a  tangential  electric  field  at  the  sur¬ 
face  of  the  antenna.  The  electric  field,  together  with  the  known  magnetic 
field  at  the  wire  surface,  represents  a  flow  of  complex  power  normal  to 
the  wire  surface.  This  complex  power  is  used  to  define  the  antenna  im¬ 
pedance.  For  the  true  current  distribution  there  is,  of  course,  no  flow 
of  power  normal  to  the  wire  surface  except  at  the  gap  where  the  emf 
exciting  the  antenna  is  inserted.  The  accuracy  of  the  Poynting-vector 
method  and  of  the  emf-method  depends  upon  the  accuracy  of  the  assumed 
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current  distribution.  As  one  might  expect,  the  accuracy  is  usually  much 
poorer  for  the  reactance  than  for  the  resistance  because  the  reactance 
depends  much  more  strongly  on  the  near-zone  field  than  does  the  resis¬ 
tance,  The  accuracy  usually  is  reasonably  good  for  thin,  resonant 
dipoles  for  which  the  current  distribution  can  be  predicted  with  good 
accuracy.  The  accuracy  5.s  not  very  good  for  general  wire  antennas, 
because  the  current  distribution  cannot  be  predicted  with  the  required 
precision.  The  emf-method  was  originally  developed  by  Brillouin^ 
in  1922. 

Both  the  Poynting -vector  method  and  the  emf-  method  make  use 
of  an  assumed  current  distribution  in  the  calculation  of  the  antenna 
impedance.  The  first  attempt  to  obtain  a  m^re  accurate  expression 
for  the  impedance  by  actually  solving  for  the  current  distribution  on  a 
wire  antenna  was  made  by  Hallen^  who  derived  an  integral  equation  for 
the  current  on  a  thin  cylindrical  dipole.  Hallen,s  analytical  solution 
of  this  integral  equation  has  been  improved  by  an  iteration  procedure 
and  gives  acceptable  results  for  a  thin  dipole. 

By  considering  the  dipole  to  be  a  spherical  transmission  line, 
Schelkunoff  5  could  express  the  total  field  around  a  dipole  as  a  sum 
of  spherical  modes.  For  a  dipole  of  arbitrary  shape,  these  modes  are 
coupled  and  the  solution  for  the  mode  amplitudes  satisfy  an  infinite  set 
of  coupled  differential  equations.  For  the  thin  dipole,  certain  approxi¬ 
mations  are  possible  because  the  principal  spherical  TEM-mode  is 
excited  much  more  strongly  than  other  modes,  and  a  reasonably  ac¬ 
curate  expression  for  the  input  impedance  can  be  obtained.  In  the 
special  case  of  a  biconical  antenna,  an  analytical  solution  can  be  ob¬ 
tained  because  the  principal  mode  coupling  occurs  over  a  limited  region 
near  the  ends  of  the  antenna. 

The  development  of  the  Schwinger  variational  principle  during 
World  War  II  has  enabled  the  approximate  solution  of  numerous  com¬ 
plicated  electromagnetic  boundary- value  problems.  The  cylindrical 
dipole  is  among  the  problems  of  practical  interest  that  have  been  treated 
by  this  method.  Using  Hallen’s  integral  equation  for  the  current  dis¬ 
tribution  on  a  cylindrical  dipole,  Tai ^  derived  a  variational  expression 
for  the  input  impedance.  Tai  showed  that  the  assumption  of  a  simple 
standing  wave  of  current  on  the  dipole  leads  to  the  same  impedance  as 
the  emf  method.  More  accurate  results  could  be  obtained  by  using  more 
complex  current  distributions, 

A  review  of  the  last  eighty  years  of  progress  in  analytical  work 
on  the  linear  antenna  has  been  given  by  R,  W.  King,^ 

The  limitations  of  the  analytical  method  described  above  can  be 
overcome  by  solving  the  integral  equations  for  the  antenna  current 
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distribution  directly  by  numerical  techniques.  This  approach  has  been 
used  by  Mei®  to  solve  a  voltage-type  integral  equation  (derivable  from 
a  voltage,  analogous  to  Hallen's  integral  equation)  for  the  current  dis¬ 
tribution  on  a  wire  antenna  of  arbitrary  geometry.  Mei  presented 
numerical  results  for  the  current  distribution  on  dipole,  loop,  and 
spiral  antennas.  Similar,  but  less  general,  techniques  have  been  used 
by  Baghdasarian  and  Angelakos  9  to  calculate  the  current  distribution 
on  loop  antennas,  and  by  Richmond  10  to  calculate  the  scattering  from 
a  straight  wire.  Most  recently,  Harrington  H  has  solved  an  electric- 
type  integral  equation  for  the  current  distribution  on  a  wire  antenna. 
Harrington  presented  numerical  results  for  a  dipole. 

Neither  the  voltage-type  integral  equation  nor  the  electric- 
type  integral  equation  are  ideally  suited  for  numerical  calculations. 

The  voltage-type  integral  equation  does  treat  sources  and  sinks  along 
the  antenna  with  good  accuracy,  but,  because  it  invdves  a  double  inte¬ 
gration  over  the  antenna,  it  is  more  complicated  than  the  electric-type 
integral  equation.  On  the  other  hand,  the  electric-type  integral  equation 
from  which  the  current  distribution  can  be  found  subject  to  conditions  on 
the  electric  field  along  the  wire  surface,  is  not  very  convenient  to  use 
near  sources  and  sinks  where  the  tangential  electric  field  will  vary  quite 
rapidly.  In  the  present  report,  these  difficulties  have  been  overcome  by 
using  a  combination  of  an  electric-type  and  a  voltage-type  integral  equa¬ 
tion,  The  electric-type  integral  equation  is  used  to  satisfy  the  boundary 
conditions  for  the  electric  field  at  the  surface  of  the  wire  composing  the 
antenna  (excluding  the  gaps).  At  the  gaps  of  the  antenna  where  sources 
or  sinks  may  be  located,  a  voltage-type  equation  is  satisfied. 

As  shown  in  Fig.  1,  the  wire  antenna  considered  in  this  report 
may  contain  an  arbitrary  number  of  gaps,  all  of  whicn  can  be  coupled 
together  by  an  arbitrary  non-radiating  network;  e.g.  ,  by  transmission 
lines  as  they  are  used  in  log-periodic  antennas.  When  the  antenna  gaps 
are  coupled  by  a  network,  the  gap  conditions  also  become  coupled 
equations  in  the  gap  currents. 
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Arbitrary  Wire  Antenna 
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METHODS  OF  CHECKING  THE  ACCURACY 


When  one  develops  a  computer  program  to  calculate  the  proper¬ 
ties  of  virtually  any  wire  antenna  with  very  high  accuracy,  it  is  im¬ 
portant  to  be  able  to  check  out  the  program  by  comparison  with  results 
that  have  been  obtained  by  other  methods.  The  only  theoretical  case 
for  which  a  comparison  is  worth  while  is  the  simple  dipole  antenna. 

In  most  cases  one  must,  therefore,  rely  on  experimental  measurements 
for  a  comparison.  The  quantity  best  suited  for  a  comparison  is  the  in¬ 
put  impedance  of  the  antenna.  The  input  impedance  describes  an  impor¬ 
tant  property  of  the  antenna  in  terms  of  only  two  numbers  --  a  resistance 
and  a  reactance.  From  an  experimental  point  of  view,  the  input  im¬ 
pedance  can  be  determined  with  much  higher  accuracy  than  either  the 
current  distribution  or  the  radiation  pattern  of  the  antenna.  Also,  from 
a  theoretical  point  of  view,  it  is  known  that  the  input  impedance  calcu¬ 
lated  at  the  gaps  of  a  wire  antenna  will  often  depend  quite  critically 
upon  the  accuracy  of  the  current  distribution  while  the  radiation  pattern 
is  usually  less  sensitive  to  small  errors  in  the  current  distribution. 

The  thin  center-fed  half-wave  dipole  is  an  excellent  example  of  this 
fact.  The  current  distribution  most  often,  used  on  this  dipole  is  of 
sinusoidal  form  --  a  distribution  which  is  extremely  close  to  the  cor¬ 
rect  form  and  which  will  also  permit  the  radiation  pattern  to  be  cal¬ 
culated  with  high  accuracy.  However,  the  input  impedance,  which  is 
proportional  to  the  derivative  of  the  current  at  the  gap,  is  zero  when 
the  current  distribution  is  assumed  to  bo  sinusoidal.  The  current  dis¬ 
tribution  and  radiation  pattern  will,  therefore,  often  be  calculated  with 
much  higher  accuracy  than  the  input  impedance. 

From  the  foregoing  discussion,  we  conclude  that  the  input  im¬ 
pedance  will  often  be  better  suited  for  checking  calculations  than  either 
the  current  distribution  or  the  radiation  pattern.  This  may  not  be  true 
in  cases  where  parts  of  the  antenna  are  very  weakly  coupled  to  the  input 
terminals  of  the  antenna  and  where,  therefore,  the  current  distribution 
could  conceivably  possess  errors  that  would  not  show  up  very  strongly 
in  the  input  impedance.  Such  errors,  if  present,  would  of  course  affect 
the  accuracy  of  the  radiation  pattern.  Whether  such  errors  are  present 
or  not  can  be  detei mined  by  calculating  the  radiation  resistance  of  the 
antenna  by  the  Poynting-vector  method:  integrating  the  radiation  pattern 
over  all  directions  of  space.  If  the  radiation  resistance  thus  obtained 
does  not  agree  with  the  resistance  calculated  at  the  input  terminals  of 
the  antenna,  the  errors  mentioned  are  probably  present.  Most  of  the 
numerical  results  presented  in  this  paper  are  input  impedance  data. 

When  the  antenna  is  considered  as  a  transmitting  antenna,  the 
source  of  the  current  distribution  induced  in  a  wire  antenna  is,  of  course, 
emfTs  introduced  at  various  places  in  the  antenna.  Quite  often,  however, 
the  antenna  is  considered  a  receiving  antenna,  and  the  problem  is  now 
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to  calculate  the  current  distribution  induced  in  the  wire  antenna  by  an 
incident  plane  wave,  and  to  calculate  the  scattering  pattern  of  the 
receiving  antenna.  The  source  is,  in  this  case,  the  incident  plane 
wave.  The  computer  program  developed  for  the  general  wire  antenna 
accepts  both  types  of  sources. 


3. 


THE  INTEGRAL  EQUATIONS  FOR  ARBITRARY  WIRE  ANTENNAS 


The  electromagnetic  field  radiated  by  a  wire  antenna  is  given  in 
terms  of  an  integral  transform  of  the  current  distribution.  Therefore, 
the  evaluation  of  the  radiation  field  is,  in  general,  a  simple  matter 
once  the  current  distribution  is  known.  The  current  distribution  on  the 
antenna  is  correct  when  the  electric  field  it  radiates  has  zero  tangential 
components  on  all  metallic  surfaces  on  the  antenna  and  when  the  current 
distribution  in  the  gaps  satisfies  certain  conditions  to  account  for  the 
sources  and  sinks  placed  there.  Satisfying  these  boundary  conditions 
at  the  antenna  surface  and  in  the  gaps  leads  to  an  integral  equation  for 
the  current  distribution,  which  we  shall  now  derive. 

The  time  factor  e  Wt  and  m.k.s.  units  will  be  used  throughout 
the  report. 

The  electric  field  vector  E  produced  by  a  current  density  dis¬ 
tribution  J  in  free  space  can  be  expressed  as  a  volume  integral  over 
the  source  distribution 


E(r»)*-  J(r)*G(?,r*)dv. 

V 


(1) 


Here,  r*  is  the  radius  vector  of  the  observation _point,  and  r  is  the 
radius  vector  of  the  source  point.  The  function  G  is  a  dyadic  which, 
according  to  Silver,  *2  is  given  by 


G  (r,  r1)  * 


fx(kR}RR+  f2(kR)  e 


,  -jkR 


kR 


(2) 


where  k  *  2tt/X  is  the  free  space  propagation  constant.  R  is  the  dis¬ 
tance  from  the  source  point  to  the  observation  point,  R  is  a  unit  vector 
pointing  from  the  observation  point  to  the  source  point  (see  Fig.  2),  and 
I  is  the  unit  dyad.  The  functions  f^  and  f^  are  defined  as  follows: 


(x)  * 

»l  +  j  -+  -2  . 

X  X 

(3) 

(x)  * 

1  . 1  1 

1  -  J  X  ’  x2  * 

(4) 
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Let  us  apply  Eq.  (1)  to  a  wire  antenna.  In  a  wire  antenna,  the 
current  flows  on  the  surface  of  wires.  For  simplicity,  we  shall  assume 
in  the  following  that  the  current  flows  on  the  axis  of  the  wire.  The  rela¬ 
tive  error  thereby  introduced  in  calculations  of  the  electric  field  on  the 
wire  surface  is  of  the  order  of 

1  -  J^kr^)  ~  (kr^)2  where  r^  is  the  wire  radius. 

For  most  practical  purposes,  this  error  is  negligible.  Correction  terms 
compensating  for  this  error  could  be  derived,  if  necessary. 

Assuming  the  current  vector  on  the  axis  of  the  wire  to  be  I(r), 
Eq.  (1)  reduces  to  the  following  line  integral  over  the  axes  of  the  wire 
structure 


E(P) 


{t}*  G(r,  r‘)d(kt) 


C 


(5) 


where  t  is  a  coordinate  measured  along  the  wire,  I(t)  is  the  current 
flowing,  and  i(t)  is  the  unit  vector  in  the  defined  positive  direction 
of  current  flow  at  the  coordinate  t  (see  Fig.  2),  The  integration  is 
extended  over  the  entire  wire  structure.  When  Eq.  (2)  for  G  is 
introduced  in  this  equation,  we  obtain 


E(P) 


l(t)  f1(kR)i(t)'RR+f2  (kR)i(t) 


5“jkR 

kR 


d  (kt),  (6) 


To  arrive  at  the  integral  equation  for  the  current  distribution  I  (t), 
we  shall  now  satisfy  the  boundary  condition  for  the  tangential  electric 
field  at  the  surface  of  the  wire.  We  shall  assume  that  the  source  of  the 
current  distribution  induced  on  the  wire  structure  includes  an  external 
source  which  radiates  a  wave  incident  on  the  antenna  with  the  electric 
field  vector  Einc(t*)  at  coordinate  t*  of  the  antenna.  The  boundary 
condition  at  the  surface  of  the  wire  requires  that  the  total  tangential 
electric  field  equal  the  product  of  the  internal  wire  impedance  per  unit 
length  Zw  and  the  current  flowing  in  the  wire.  The  total  field  at  the 
wire  surface  is  the  sum  of  the  incident  field  and  the  field  expressed  by 
Eq,  (6).  We  find,  therefore,  the  following  integral  equation  for  the 
current  distribution  I  (t) 
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The  observation  point  at  which  the  boundary  condition  is  satisfied 
has  the  coordinate  t*  measured  along  the  wire  structure  and  is  located 
a  distance  from  the  wire  axis  equal  to  the  wire  radius  at  the  point  t*, 
while  t*  is  a  unit  vector  pointing  along  the  wire  surface  at  the  obser  vation 
point  (t*  has  no  circumferential  component  to  the  wire),  as  shown  in  Fig.  2. 
We  assume  here  that  the  wires  composing  the  antenna  are  of  conical  shape, 
t*  is  therefore  not  necessarily  identical  to  i  (t*)  but  does,  in  general, 
point  in  a  different  direction  depending  upon  where  the  observation  point 
is  chosen  along  the  wire  perimeter.  Strictly,  Eq.  (7)  should  hold  for  any 
point  on  the  surface  of  the  wire.  However,  the  equation  cannot  be  satis¬ 
fied  in  all  points  along  the  circumference  of  the  wire  because  we  have 
tacitly  neglected  circumferential  currents  on  the  wires,  an  approximation 
which  is  extremely  good  in  most  practical  cases  and  is,  of  course,  better 
the  thinner  the  wire  is.  Accordingly,  it  should  make  little  difference 
where  we  choose  the  observation  point  on  the  wire  circumference.  Yet, 
having  neglected  the  circumferential  current  component,  we  can  some¬ 
times  do  better  than  just  choosing  an  arbitrary  observation  point  on  the 
wire  circumference.  Considering,  for  example,  the  incident  field  on  the 
left-hand  side  of  Eq.  (7),  it  is  more  accurate  to  evaluate  the  field  in  the 
direction  of  the  wire  axis  in  an  observation  point  placed  on  the  wire  axis 
itself  than  it  is  to  evaluate  the  field  in  any  point  on  the  wire  circumference. 
On  the  left-hand  side  of  Eq.  (7),  we  therefore  replace  t1  by  i(t’)  and 
tacitly  assume  the  observation  point  to  be  on  the  wire  axis.  On  the  right- 
hand  side  of  Eq,  (7),  we  shall  also  replace  t*  by  i  (t*)  and  assume  the 
observation  point  to  be  on  the  wire  axis  rather  than  on  the  wire  surface. 

We  require,  however,  that  the  wire  axis  at  the  source  point  and  the  wire 
axis  at  the  observation  point  do  not  coincide.  If  the  wire  axes  do  coincide, 
the  observation  point  is  chosen  on  the  wire  surface.  Since  the  field 
radiated  by  the  source  is  in  that  case  symmetric  about  the  wire  axis,  the 
observation  point  can  be  chosen  anywhere  on  the  wire  circumference. 
Accordingly,  on  the  right-hand  side  of  Eq.  (7),  let  us  replace  t*  by  a 
vector  s'  defined  as  follows 


i  (t*) 


when  R  x  i  (t1)  ^  0 
when  Rxi  (t*)  *  0. 


(8) 
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Tacitly  assuming  the  observation  point  to  be  chosen  as  described  above, 
after  a  few  transformations  of  Eq.  (7)  we  arrive  at  the  following  integral 
equation  for  I  (t) 


i?(t*).EinC(t*)»  ~^I(t*) 


j  3 0  s^I  (t)  jf  l (kR)  i  (t) *  RR  +  f2(kR)  i  (t) 


;-jkR 

kR 


d(kt). 


(9) 
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4. 


REDUCTION  TO  A  SET  OF  LINEAR  EQUATIONS 


In  order  to  obtain  a  solution  of  the  integral  equation,  Eq,  (9),  we 
shall  reduce  it  to  a  set  of  simultaneous  linear  equations  in  sampled 
values  of  the  current  distribution.  Assume  the  total  length  of  wire  to 
be  divided  into  N  intervals.  Let  the  sampled  value  of  the  current  in 
the  center  point  of  the  n*th  interval  ^  be  1^  and  let  the  coordinate 
in  that  point  be  tn  (see  Fig.  1).  The  integral  in  Eq.  (9)  can  now  be 
approximated  by  a  linear  combination  of  the  N  sampled  currents,  and 
these  currents  can  be  solved  for  by  satisfying  the  integral  equation  (9) 
in  the  N  sampling  points.  This  leads  to  a  set  of  N  linear  equations 
which  can  be  solved  by  a  digital  computer.  The  accuracy  of  the  solution 
will  depend  quite  strongly  upon  how  the  reduction  of  the  integral  in  Eq, 

(9)  to  a  linear  combination  of  the  sampled  current  values  is  accomplished. 
It  is  important  that  the  final  current  distribution,  calculated  to  satisfy  the 
boundary  conditions  correctly  in  the  sampling  points  chosen,  also  pro¬ 
duce  a  boundary  field  that  varies  smoothly  between  the  sampling  points. 
Otherwise,  the  boundary  field  which  acts  as  a  virtual  emf  that  drives  an 
erroneous  current  through  the  antenna,  could  significantly  impair  the 
accuracy  of  the  solution.  We  have  found  that  a  constant  or  linearly 
varying  current  in  each  interval  of  the  antenna,  as  used  by  Harrington,  ^ 
can  lead  to  significant  errors  when  the  interval  length  is  much  greater 
than  the  wire  diameter.  Instead,  the  following  representation  for  the 
current  distribution  in  the  n*th  !nterval  of  the  antenna  was  found  to  lead 
to  very  satisfactory  accuracy  in  most  cases: 


I (t)  =  T  +  C  cos  k(t-t  }  +  S  sin  k(t~t  )  .  (10) 

w  n  n  n  n  n 


This  distribution  represents  a  sum  of  a  constant  current  T^  and  a  stand¬ 
ing  wave  with  propagation  constant  of  free  space.  This  distribution  has 
been  used  also  by  Mei.  ® 

When  Eq.  (10)  is  introduced  in  the  integral  equation  Eq.  (9),  and 
Eq,  (9)  is  required  to  be  satisfied  in  the  N  sampling  points,  we  find, 
by  integrating  over  the  N  integration  intervals  the  following  set  of 
equations : 
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where  x  ,  E  ,  Z  ,  8  ,  and  I  are  the  values  of  i,  E  ,  Z  ,  s, 

and  I  at^e  m%i  sampling  point,  ancP  An  is  the  length  of  the  n*th  inte¬ 
gration  interval,  Rnm  is  the  vector  pointing  from  the  n*th  sampling 
point  to  the  m*th  observation  point,  and  in  is  the  value  of  i  at  the  n’th 
sampling  point. 


The  vector  functions  G0,  GQ,  and  Gs  appearing  in  Eq.  (11) 
represent  the  electric  field  radiated  by  a  dipole  of  length  An,  oriented 
in  the  direction  in,  and  carrying  a  constant,  cosine-  or  sine-distribu¬ 
tion  of  current  respectively.  The  field  is  evaluated  in  a  point  with  radius 
vector  R_m  pointing  from  the  center  of  the  dipole.  The  functions  GQ, 
Gc,  and  Crs  are  derived  in  Appendix  A, 


As  shown  in  Appendix  B,  the  constants  Tn,  Cn,  and  Sn  can  be 
expressed  as  a  linear  combination  of  a  maximum  of  three  sampled  values 
of  current  because  a  constant  -  cosine  -  sine  curve  can  be  fitted  to  three 
points  only.  When  ^n*  ^n'  and  Sn  can  be  expressed  as  linear  com¬ 

binations  of  the  current  values  in  the  N  sampling  points,  Eq.  (11)  ob¬ 
viously  can  be  interpreted  as  a  set  of  N  linear  equations  which  can  be 
solved  for  the  N  unknown  sampled  current  values.  If  we  consider  the 
antenna  as  a  transmitting  antenna,  however,  the  incident  field  Einc  is 
zero,  and  the  solution  of  Eq.  (11)  is  zero,  eigensolutions  of  the  equation 
being  ignored.  The  reason  why  the  solution  turns  out  to  be  zero  is  ob¬ 
viously  that  we  have  not  yet  defined  how  we  want  to  excite  the  antenna. 

An  additional  condition  must  be  satisfied  at  each  gap  of  the  antenna.  This 
gap  condition  can  also  be  reduced  to  a  linear  combination  of  the  sampled 
current  values.  However,  for  each  gap  condition  added  to  Eq.  (11),  we 
have  to  introduce  an  extra  unknown  to  keep  the  total  number  of  unknowns 
the  same  as  the  total  number  of  equations.  The  most  natural  choice  for 
this  unknown  is  the  gap  current  itself.  When  gaps  are  present,  some  of 
the  constants  Tn,  Cn,  and  Sn  will  also  depend  upon  the  gap  currents. 
Assuming  the  presence  of  Nga  5  gaps,  and  numbering  the  gap  currents 
from  N+  1  through  N+Ngap,  jfcq,  (11)  can  therefore  be  written  in  the 
following  abbreviated  form: 
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(12) 


The  gap  conditions  are  quite  complex  in  the  general  case  where 
several  wires  join  at  the  gap,  where  sources  and/or  sinks  can  be  in¬ 
serted  into  each  wire  of  the  junctions,  and  where  all  junctions  can  be 
connected  by  an  arbitrary  non-radiating  netwoi'k,  as  shown  in  Fig.  1. 

In  order  to  understand  the  nature  of  the  gap  condition,  we  shall  set  up 
the  condition  in  a  simplified  situation.  Assume  that  each  gap  joins 
only  two  wires,  and  that  these  gaps  are  not  connected  by  a  network. 
Assume  furthermore  that  the  wires  are  conical  near  the  gap  (see  Fig.  3), 
The  dominant  field  near  the  gap  can  then  be  represented  by  a  TEM-wave 
on  a  biconical  transmission  line,  and  the  standard  transmission  line 
equations  permit  establishing  the  desired  relationship  between  the 
current  distribution  and  the  source  generating  it.  Accordingly,  we 
find  for  the  terminal  voltage  at  the  gap 
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d(kt) 
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(13) 


where  is  the  generator  emf,  Z^  is  the  generator  impedance,  ZQ 
is  the  characteristic  impedance  of  the  biconical  transmission  line,  and 

dl 

d(kt) 

is  the  derivative  of  the  current  on  the"positive"  side  of  the  gap  (t >  ) . 

Equation  (13)  is  correct  only  if  the  current  distribution  is  symmetric 
about  the  gap.  This  cannot,  in  general,  be  assumed  to  be  the  case. 

The  current  distribution  will,  therefore,  in  general  also  contain  a  first- 
order  spherical  mode  near  the  gap.  While  this  mode  does  not  contribute 
to  the  gap  current  itself,  it  does  contribute  to  the  derivative  of  the  cur¬ 
rent.  Fortunately,  however,  the  derivative  of  the  TEM-wave  is  an  odd 
function  about  the  gap,  while  the  first-order  wave  is  an  even  function 
about  the  gap.  The  first-order  wave  can,  therefore,  be  eliminated  from 
Eq.  (13)  by  replacing  the  derivative  in  Eq.  (13)  by  half  the  difference 
between  the  derivatives  on  the  two  sides  of  the  gap 
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The  derivates  in  this  equation  can,  of  course,  be  approximated 
numerically  quite  well  as  a  linear  combination  of  a  few  of  the  sampled 
currents,  including  the  gap  current.  Omitting  the  details  of  this  cal¬ 
culation,  \va  find  that  for  Ngap  gaps,  Eq.  (14)  can  be  reduced  to  the 
following  set  of  equations: 


N+N 


I 

n  *  1 


gap 

H  I  *  V.  ;  m  a  N+  1,  N+ 2, 
mn  n  1m 


•»N+N  .  (15) 

gap 


where  Vjm  is  the  emf  in  the  gap  numbered  m. 


Equations  (12)  and  (15)  form  a  set  of  N+Ngap  simultaneous  linear 
equations  from  which  the  N+  N_a_  unknown  current  values  can  be  found. 
These  equations  can  be  solved  Dy  a  digital  computer.  When  the  currents 
have  been  found,  the  input  impedance  of  the  antenna  seen  from  the  various 
gaps  is,  of  course,  immediately  found  as  the  ratio  between  the  voltage 
across  the  gap  and  the  gap  current 


Zinput  *  Vi5?-"3”i52.  .  m  at  N+  1,  N+2 _ _  N+N  (16) 

m  I  gap 

m 


where  Zjm  is  the  generator  impedance  at  the  gap  numbered  m  and 

I  is  the  gap  current, 
m 
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THE  RADIATION  PATTERN 


When  the  unknown  currents  have  been  found,  it  is,  of  course, 
also  an  easy  matter  to  calculate  the  radiation  pattern.  The  radiated 
far-zone  field  can  be  calculated  basically  from  the  same  set  of  for¬ 
mulas  derived  in  Appendix  A,  but  these  formulas  can  be  simplified 
considerably  when  the  observation  point  is  in  the  far  zone.  Without 
going  into  the  details  of  these  calculations,  we  find  the  following  ex¬ 
pression  for  the  electric  field  vector  E  in  the  far  zone 


i  E  *  j  30 


-jkR 
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E  +  C  E  +  S  E 
on  n  cn  n  sn 


where  R  is  the  distance  from  the  origin  of  the  coordinate  system  to 
the  observation  point  in  the  far-zone,  and  R  is  a  unit  vector  pointing 
in  that  direction.  Rn  is  Jhe  radius  vector_of  the  center  of  the  n*th 
interval  of  the  antenna.  EQn,  Ecn,  and  Esn  have  components  in  the 
direction  of  i  and  of  the  perpendicular  unit  vector 


which  is  normal  to  i  , 

n 


For  these  components  the  following  expressions  are  found: 
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(22) 


These  expressions  are  easily  evaluated  by  a  digital  computer. 
Since  the  input  resistance  of  the  antenna  is  already  known,  the  absolute 
gain  of  the  antenna  can  be  calculated.  The  far-zone  field  will  generally 
be  elliptically  polarized.  The  computer  program  developed  calculates 
the  gain  along  the  two  principal  axes  of  the  polarization  ellipse,  the 
total  gain,  the  axial  ratio,  the  orientation,  and  the  sense  of  rotation  of 
the  polarization  ellipse. 
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SYMMETRY  CONSIDERATIONS 


Often,  an  antenna  will  possess  certain  symmetry  properties 
such  that  only  part  of  the  total  current  distribution  is  unknown  while 
the  rest  of  the  current  distribution  can  be  described  in  terms  of  the 
unknown  current  distribution  (e.g. ,  an  antenna  and  its  image  in  a 
perfect  ground  plane}.  A  considerable  amount  of  computer  time  can 
be  saved  by  taking  these  symmetry  properties  into  account  and  solving 
for  only  that  part  of  the  current  distribution  that  is  really  unknown. 

In  setting  up  the  equations  for  the  unknown  current,  distribution,  it  is, 
of  course,  necessary  to  take  the  mutual  coupling  between  the  unknown 
current  distribution  and  the  rest  of  the  current  distribution  properly 
into  account. 

The  computer  program  developed  will  utilize  the  symmetry 
properties  of  an  antenna  and  its  image  in  a  perfectly  conducting 
ground  plane,  and  will  also  utilize  the  symmetry  properties  of  the 
azimuthal  modes  of  a  ring  array  of  an  arbitrary  number  of  equispaced 
wire  antennas. 
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7. 


NUMERICAL  RESULTS 


The  computer  program  developed  will  solve  for  the  current 
distribution,  the  input  impedance-  and  the  radiation  pattern  of  virtually 
any  ring  antenna  array  whose  basic  element  is  a  wire  antenna  of  essen¬ 
tially  arbitrary  geometry.  There  is  no  basic  limitation  to  the  electrical 
size  of  the  antenna  that  can  be  handled,  but  the  cost  of  the  calculations 
increases  at  a  rate  between  the  square  and  the  cube  of  the  electrical  size. 
The  maximum  integration  interval  length  that  can  be  chosen  depends  upon 
geometry,  but  can  seldom  be  allowed  to  exceed  0.  25  wavelengths  with¬ 
out  a  considerable  loss  of  accuracy.  When  the  matrix  for  the  equations 
to  be  solved  by  the  computer  does  not  fit  into  the  memory  of  the  com¬ 
puter,  part  of  the  matrix  will  be  on  magnetic  tape. 

To  illustrate  the  use  of  the  computer  program  developed,  we 
shall  now  show  some  numerical  results  obtained  for  a  number  of  wire 
antennas.  The  results  have,  in  most  cases,  been  verified  by  measure¬ 
ments. 


The  only  wire  antenna  that  it  has  been  possible  to  analyze  with 
high  accuracy  by  analytical  methods  is  the  simple  dipole.^*  ^  Since 
the  analytical  results  for  the  dipole  are  supposedly  more  accurate  than 
most  measured  results,  it  is  very  valuable  to  be  able  to  check  the 
accuracy  of  the  general  computer  program  developed  by  comparison 
against  the  analytical  results  for  a  dipole.  Table  1  shows  the  values 
of  input  impedance  for  a  cylindrical  dipole  obtained  by  the  second-order 
approximation  method  of  King  and  Middleton.  ^  The  Table  also  shows 
the  results  obtained  by  the  present  computer  program  when  the  dipole 
is  excited  through  a  conical  gap  extending  over  ten  percent  of  the  length 
of  the  dipole.  The  agreement  between  the  values  of  input  resistance 
is  seen  to  be  extremely  good.  Some  differences  are  observed  in  the 
reactance  values,  probably  due  to  the  fact  that  the  dipole  treated  by  the 
present  program  has  a  conical  gap.  The  total  number  of  integration 
intervals  used  for  the  calculations  in  Table  1  was  20,  In  most  practi¬ 
cal  cases  where  the  accuracy  need  not  be  nearly  so  high  as  in  Table  1, 
many  fewer  intervals  can  be  used. 

In  most  cases,  the  computer  program  developed  for  calculating 
the  properties  of  arbitrary  wire  antennas  must  be  checked  out  by  care¬ 
ful  experimental  measurements.  Figures  4  through  12  show  the  cal¬ 
culated  and  measured  input  impedance  of  various  wire  antennas. 

Figures  4  and  5  show  the  input  impedance  of  a  fat  monopole  and  of  a  fat 
conical  monopole.  The  agreement  between  the  theoretical  and  the  ex¬ 
perimental  results  is  seen  to  be  excellent. 


Table  1 


Input  Impedance  of  Linear  Antenna  above  Perfect  Ground 
Antenna  Height  h.  Wire  Radius  a. 


Resistance,  Ohms 

Reactance 

,  Ohms 

kh 

Andreas en- 
Tanner 

King- 

Middleton 

Andreas en- 

Tanner 

King- 

Middleton 

0.5 

2.46 

2.50 

-  639 

-  628 

0.7 

5. 14 

5. 14 

-  409 

-  405 

0.9 

9.07 

9.07 

-  266 

-  267 

1.1 

14.8 

14.7 

-  162 

-  165 

1.3 

22.9 

22.8 

-  75. 1 

-  81.8 

1.5 

34.9 

34.7 

3.9 

-  5.1 

1.7 

53 .0 

52,9 

82.2 

71.4 

1.9 

81.5 

82,5 

166 

155 

-■■■■— — — . — 
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Fig.  4.  Input  Impedance  of  Fat  Monopole  above  Ground  from 
80  MHz  to  120  MHz  in  Steps  of  10  MHz 
(Chart  Impedance  *  50  ohms) 
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Fig.  5.  Input  Impedance  of  Conical  Monopole  above  Ground 
from  80  MHz  to  200  MHz  in  Steps  of  20  MHz 
(Chart  Impedance  *  50  ohms) 
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Figure  6  shows  a  2.5  to  1  log-periodic  array  of  five  elements 
fed  from  a  200  ohm  transmission  line.  The  calculated  and  measured 
input  impedance  of  this  antenna  over  the  2.5  to  1  frequency  band  in 
which  it  should  be  operative  is  shown  in  Figs.  7  and  8.  The  agreement 
betweenlhe  theoretical  and  the  measured  results  is  seen  to  be  ex¬ 
tremely  good.  Even  the  smaller  loops  of  the  impedance  curve,  which 
are  due  to  the  reflection  from  the  open-ended  transmission  line,  agree 
very  well. 

In  the  above  examples,  we  have  shown  only  the  calculated  input 
impedance  data.  The  computer  output  also  includes  the  current  dis¬ 
tribution  and  the  radiation  pattern  of  the  antennas  investigated,  but  it 
would  be  impossible  to  show  all  of  this  information  in  this  report. 

Figure  9  illustrates  the  calculation  of  a  radiation  pattern  of  a  ten- 
element  Yagi  antenna  array  previously  calculated  by  King.^  The 
radiation  pattern  calculated  using  the  present  program  is  seen  to  agree 
well  with  King*s  calculated  pattern. 

With  the  wire  antenna  program  we  can  treat  with  surprisingly 
good  results  some  configurations  that  one  does  not  at  first  associate 
with  linear  antennas. 

Figure  10  shows  an  example  of  a  structure  composed  of  cylinders, 
tapered  cylinders  and  cones,  that  has  more  than  a  passing  resemblance 
to  an  aircraft  with  a  tail-cap  antenna.  In  this  example,  all  of  the 
members  —  including  the  air-foil  surfaces  --  have  circular  cross 
sections,  and  all  except  the  vertical  stabilizer  with  the  tail-cap  are 
co-planar.  The  dimensions  of  the  members  are  chosen  so  that  the 
structure  should  be  an  approximate  electrical  representation  of  the 
RB-66  aircraft.  A  much  better  representation  could  have  been  obtained 
by  using  two  or  more  parallel  cylinders  to  represent  the  flattened  air¬ 
foil  surfaces,  dropping  the  nacelles  below  the  wings  and  supporting  them 
on  struts  as  they  are  in  the  actual  aircraft  and  other  refinements.  Such 
refinements  would  increase  both  the  time  required  to  prepare  the  data 
for  the  computer  and  the  computer  time  required  to  obtain  the  desired 
solutions.  If  one  were  in  earnest  about  the  problem,  however  such 
time  would  be  well  spent.  Again,  one  could  obtain  as  much  accuracv 
as  desired  at  the  cost  of  increased  computer  time  resulting  from  the 
increased  complexity  of  the  structure  analyzed. 

In  spite  of  the  comparative  crudity  of  the  representation  shown 
in  Fig,  10,  however,  surprisingly  useful  and  accurate  results  are  ob¬ 
tained.  (The  inaccuracy,  it  should  be  pointed  out,  is  not  in  the  answers 
provided  by  the  computer  for  the  structure  analyzed,  but  in  the  repre¬ 
sentation  of  the  actual  aircraft  by  that  structure. )  The  computer  cal¬ 
culation  predicts  the  resonant  peak  in  radiation  resistance  associated 
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with  airframe  resonance  as  well  as  other  aspects  of  the  impedance  be¬ 
havior.  When  suitable  corrections  are  made  to  allow  for  the  rather 
substantial  disparity  between  the  actual  feed  gap  region  and  that  in¬ 
corporated  in  the  representation,  the  calculated  impedances  --  both 
resistance  and  reactance  --  agree  with  the  impedances  measured 
experimentally  on  a  model  of  the  actual  aircraft  within  approximately 
twenty  percent  at  those  points  checked. 

The  computer  calculation  provides  rather  interesting  data  on 
the  distribution  of  currents  on  the  airframe  members.  These  currents, 
for  a  frequency  near  airframe  resonance,  are  .hown  in  Fig.  11.  In 
the  figure,  the  airframe  members  are  represented  by  the  lines  forming 
the  axes  of  the  members  and  the  currents  are  plotted  vertically  -- 
picket-fence  fashion  --  along  these  lines.  The  phases  of  the  currents 
vary  little,  only  a  few  degrees  from  one  end  of  any  member  to  the 
other,  and  are  indicated  by  the  phases  shown  on  the  diagram,  which 
apply  to  the  currents  at  the  center  of  each  linear  member  composing 
the  airframe.  The  finite  values  of  current  at  the  extremities  of  the 
members  show  the  amount  of  curri  nt  flowing  onto  the  flat  end-caps  of 
the  cylindrical  members. 

Another  good  measure  of  the  accuracy  of  the  results  provided 
by  the  calculation  is  in  the  correspondence  of  the  measured  with  the 
calculated  radiation  patterns.  This  comparison  is  shown  in  Fig.  12 
for  a  frequency  somewhat  higher  than  the  airframe  resonance  fre¬ 
quency  (the  experimental  measurements  were  made  a  number  of  years 
ago  and  only  limited  data  were  available).  It  is  evident  that  the  com¬ 
puter  has  predicted  all  of  the  major  features  of  the  patterns.  The  only 
significant  discrepancy  is  in  the  null  below  the  aircraft  in  the  longi¬ 
tudinal  principal  plane  pattern  (top  center)  where  the  computer  predicts 
a  significantly  deeper  null  than  was  measured.  The  two  small  vestigial 
pattern  lobes  --  the  one  directed  up  and  forward,  and  the  other  rear  and 
down  on  this  pattern  --  are  predicted.  Since  the  patterns  are  changing 
quite  rapidly  with  frequency,  the  discrepancy  in  the  downward  null  may 
well  be  due  to  a  relatively  minor  misjudgment  in  the  choice  of  dimen¬ 
sions  for  the  linear  element  representation  of  the  aircraft.  It  is  almost 
certain  that  a  more  elaborate  representation  would  overcome  what  dis¬ 
crepancies  do  exist. 

Calculations  were  also  made  to  verify  measurements  made  on  a 
center-fed  whip  antenna  described  elsewhere.  16  Calculated  and 
measured  values  of  the  input  impedance  are  shown  in  Table  2,  The 
agreement  between  calculations  and  measurements  is  only  fair.  The 
discrepancy  between  the  calculated  and  the  measured  impedance  values 
is  ascribed  to  differences  in  the  experimental  and  the  theoretical  models 
used  for  the  antenna. 
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FIG.  II  CURRENT  DISTRIBUTION  ON  AIRCRAFT  WITH  TAIL-CAP  ANTENNA  AT  AIRFRAME  RESONANCE 

FREQUENCY 


COMPARISON  OF  CALCULATED  AND  MEASURED  PATTERNS 
RB-66  AIRCRAFT  WITH  TAILCAP  ANTENNA 


Calculated  and  Measured  Input  Impedance 
of  Center-Fed  Whip  Antenna 


Finally,  calculations  were  made  on  the  compact  VHF  antenna 
shown  in  Fig,  13.  The  principal  results  of  the  calculations  at  30  MHz 
are  summarized  in  Table  3.  It  is  interesting  to  notice  that  the  input 
impedance  of  the  antenna  drops  as  the'  antenna  is  raised  above  the 
ground  plane.  This  behavior  is  also  expected  from  the  theory  of 
coupled  antennas.  From  average-gain  calculations  it  is  concluded 
that  the  results  for  the  input  resistance  are  within  20  percent  of  the 
true  values.  The  small  electrical  size  of  the  antenna  puts  a  serious 
strain  on  the  accuracy  of  the  program,  the  accuracy  deteriorating  as 
the  electrical  size  of  the  antenna  approaches  zero.  The  reason  for  the 
deteriorating  accuracy  for  electrically  small  antennas  is  that  the 
constant-current  term  and  the  cosine-current  term  used  in  the  expan¬ 
sion  of  the  current  distribution  within  each  integration  interval  become 
degenerate  for  zero  length  of  the  interval.  This  quasi-degeneracy  for 
electrically  small  antennas  can  be  removed  by  using  a  different  repre¬ 
sentation  for  the  current  distribution.  The  program  was  applied  un¬ 
successfully  to  a  compact  HF  antenna  with  an  overall  size  of  0,01 
wavelengths;  the  round-off  errors  due  to  the  quasi-degeneracy  of  the 
constant-  and  the  cosine-current  terms  were  found  to  be  too  large. 


Table  3 


Calculated  Data  for  Compact  VH.b -Antenna  at  30  MHz 


Height  H  Above 
Ground  Plane 
Inches 

Input  Re¬ 
sistance, 
Ohms 

Input  Re¬ 
actance, 
Ohms 

Maximum 

Gain* 

dB 

*k 

Gain  at  Ele¬ 
vation  '3  0 
dB 

Gair? at  Ele¬ 
vation  60 
dB 

0 

0.87 

-  571 

5.30 

4.04 

-  0.77 

36 

0.61 

-  938 

6.00 

4.16 

-  1.87 

120 

0.42 

-  938 

7.61 

0.50 

-11.4 

* 


Gain  is  referenced  to  an  isotropic  radiator. 
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8. 


CONCLUSION 


The  computer  program  developed  will  evaluate  the  properties 
of  virtually  any  wire  antenna.  The  results  shown  here  have  served 
only  to  illustrate  the  applicability  of  the  program.  Besides  the  ad¬ 
vantages  of  speed,  economy,  and  flexibility  relative  to  experimental 
methods,  the  numerical  approach  provides  answers  with  a  precision 
and  a  wealth  of  detail  completely  beyond  the  reach  of  experimental 
techniques  in  most  cases.  This  detail  and  precision  can  provide  the 
basis  for  enhanced  insight  into  the  operation  of  the  antennas  analyzed, 
leading  in  most  cases  to  more  efficient,  more  economical,  and  higher 
performance  designs. 
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APPENDIX  A 


Dipole  Radiation  Fields 


The  reduced  integral  equations  for  the  wire  antenna,  Eq.  (11), 
is  given  in  terms  of  vector  functions  G0,  Gc,  and  Gg  that  repre¬ 
sent  the  electric  field  radiated  by  a  dipole  carrying  a  constant,  a 
cosinoidal,  or  a  sinusoidal  current  respectively.  These  vector 
functions  will  be  derived  here. 


Calculation  of  G  (R  ,  z,  A) 
§  o 


Consider  the  dipole  shown  in  Fig.  A-l  of  length  A  and  oriented 
in  the  direction  z.  We  want  to  calculate  the  electric  field  that  this 
dipole  radiates  in  a  point  with  radius  vector  R0  using  the  center  of 
the  dipole  as  the  origin,  when  the  current  distribution  on  the  dipole  is 
of  the  form: 

I(z)  *  C?S(kz)  (A-l) 

sin 

where  k  is  the  free  space  propagation  constant.  This  problem  was 
solved  almost  a  half  century  ago  by  Brillouin.  ^  His  results  can  be 
found  in  a  convenient  form  in  the  textbook  by  Stratton.  ^  Thus,  for  the 
dipole  considered  here,  we  find  the  following  exact  expressions  for  the 
z-  and  the  p-components  of  the  electric  field 
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where  integrations  are  extended  over  the  length  of  the  dipole.  i)j  is  the 
spherical  wave  function 
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(A-7) 


The  validity  of  these  expressions  can  be  verified  very  easily.  Let  us 
evaluate  first  the  p- component  of  the  field.  Introducing 
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In  eq,  (A—  5),  we  find 
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For  the  z-component  of  Gq  we  find  from  Eq.  (A-6) 
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The  integral  over  i)/  in  this  equation  cannot  be  evaluated  exactly  in 
closed  form.  However,  a  reasonably  accurate  expression  valid  for 
small  values  of  kp  can  be  obtained  as  follows.  Let  us  divide  the 
integration  path  from 


kA 
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to 


path  from  kz  to 


kA 
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kA 
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into  one  path  from 
We  find  then 
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where 


Sign  (x) 


and  where  we  have  made  use  of  the  following  relation 


dR  t  -  s,  .  Vr2  -  p2 
dt  R  R 


(A- 14)  I 

s 


(A-i5) 


The  upper  sign  here  applies  for  t  >  z;  the  lower  sign  for  t  <  z. 


Eq.  (A- 13)  can  be  further  reduced  by  introducing  the  substi- 


whereby  we  find 


+  kA 
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where  F  (kR) 


-jkpx 


x  -  1 


(A- 17) 


The  evaluation  of  the  function  F  is  different  depending  on  whether 

kR  .  ,  kR  .  , 

is  large  or  not.  If  — —  is  large,  the  function  F  is  most  easily 

evaluated  by  expressing  it  as  follows: 


c?\  r  "jkR 

F(kR)  *  * j  2  H°  tM-  v5#^rd(kR) 


(A- 18) 


where  we  have  made  use  of  the  following  integral  representation  of  the 
Hankel  function  of  the  second  kind  and  of  order  zero: 
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The  integral  in  Eq.  (A- 18)  can  now  be  evaluated  by  expansion  of  the 
square  root  in  the  integrand  and  integrating  term  by  term.  When  we 
neglect  all  terms  of  orders  higher  than  4  in  k  p,  we  find 


F(kR)  *  -  j  |  Ho(2)(kp}+  Ci(kR)  +  j  |  -  Si(kR) 
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where  Ci  and  Si  are  the  cosine-  and  sine-integrals. 
kR 

If  —  is  not  very  large,  the  integral  over  ijj  in  Eq.  (A- 12)  is 
kp 

most  easily  evaluated  numerically;  e.g.,  by  Simpson*s  rule. 
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APPENDIX  B 


The  Current  Distribution  in  an  Integration  Interval 


In  the  reduction  of  the  integral  equations  for  the  current  distri¬ 
bution  on  the  wire  antenna  into  a  set  of  simultaneous  linear  equations 
in  sampled  values  of  the  current  distribution,  it  was  found  that  the 
accuracy  of  the  solution  depends  strongly  upon  how  the  current  distri¬ 
bution  is  assumed  to  vary  between  the  sampling  points.  The  current 
distribution  in  the  n*th  integration  interval,  as  expressed  in  Eq,  (10), 
was  found  to  yield  a  very  satisfactory  accuracy  in  most  cases.  The 
constants  Tn,  Cn,  and  Sn  can  be  expressed  as  a  linear  combination 
of  a  maximum  of  three  sampled  values  of  current,  because  a  constant- 
cosine-sine  curve  can  be  fitted  to  three  points  only.  To  find  the  con¬ 
stants  of  this  linear  combination,  assume  that  the  sampled  value  of 
current  in  the  n*th  sampling  point  with  wire  coordinate  tatn  is  In, 
and  assume  current  values  Ia  and  ib  in  neighboring  points  t  *  t 

Ct  Lf  3, 

and  t  *  t  . 

b 


If  the  curve  represented  by  Eq.  (10)  is  to  fit  the  three  points 
ta,  tr  .nd  tb,  the  constants  Tn,  Cn,  and  Sn  must  satisfy  the 
following  set  of  three  linear  equations 
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which  has  the  solution 
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where  »  sinj^k  (t^  -  t^)J  -  sin^k  (t^  -  t^)J  +  sinj^k  (t^  -  t^)J.  (B-5) 


The  currents  Ia  and  1^  in  the  above  expressions  will,  of  course, 
be  related  linearly  to  the  sampled  currents.  If  the  n*th  interval  con¬ 
sidered  above  is  an  "internal"  interval  with  at  least  one  integration  inter¬ 
val  to  each  side  of  it,  it  is  most  convenient  to  let  "a"  be  the  (n  -  l)*st 
sampling  point  and  to  let  "b"  be  the  (n  +  l)Ist  sampling  point,  thus  choosing 
Ia  *  In-i  and  1^  *  ln^  (see  Table  B-l).  There  are  a  number  of  cases, 
however,  that  require  special  consideration.  These  include  the  case 
where  one  end  of  the  n*th  interval  is  open,  and  the  case  where  the  inter¬ 
val  considered  adjoins  an  electric  or  a  magnetic  symmetry  plane.  In 
these  cases,  we  make  use  of  the  fact  that  the  current  or  its  derivative  is 
zero  at  one  end  of  the  interval  and  select  values  of  ta,  t^,  Ia,  and  1^ 
as  shown  in  Table  B-l.  Other  special  cases  listed  in  Table  B-l  include 
the  case  where  the  interval  adjoins  a  junction  of  several  wires,  and  the 
case  where  it  adjoins  a  gap  with  a  source  or  a  sink.  In  these  cases,  we 
assume  the  current  at  the  end  adjoining  the  junction  or  the  gap  to  be  un¬ 
known  (  IG  in  Table  B-l),  and  we  solve  for  this  current  by  satisfying  a 
gap  condition  as  described  in  Sec.  4. 
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APPENDIX  C 


INSTRUCTIONS  FOR  USE  OF  TRG  PROGRAM 

11  WIR  TNA  11 
Control  Data  3800  Version 


This  program  will  solve  a  very  general  class  of  problems  in  the 
area  of  radiation  from  wire  antennas  of  essentially  arbitrary  geometry. 
The  program  requires  as  input  the  coordinates  describing  the  wire 
structure,  and  information  on  sources  and  sinks  located  at  gaps  along 
the  wire  antenna,  as  well  as  information  on  non- radiating  networks  con¬ 
necting  the  gaps.  The  output  consists  of  the  current  distribution  induced 
on  the  wire  antenna,  the  antenna  input  impedance  seen  from  the  generator 
terminals  of  the  antenna,  and  the  absolute  gain  of  the  radiation  pattern  in 
any  number  of  specified  directions. 

The  program  will  handle  two  types  of  sources:  emf’s  located  in 
discrete  points  along  the  wire  structure,  and  dipole  sources  placed  any¬ 
where  with  respect  to  the  wire  structure.  The  emf-source  is  most  use¬ 
ful  in  antenna  problems,  while  the  dipole  source  makes  possible  the  solu¬ 
tion  of  problems  of  scattering  from  wire  structures. 
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Introduction 

The  program  "WIRTNA"  will  solve  for  the  current  distribution 
induced  on  wire  structures  composed  of  linear  wire  sections  joined 
arbitrarily.  The  wire  structure  is  assumed  to  be  excited  by  an  emf 
located  in  any  number  of  gaps  on  the  wire  structure,  or  by  dipole 
current  sources  placed  anywhere  with  respect  to  the  wire  structure. 

Each  linear  wire  section  is  defined  by  the  coordinates  and  by  the  wire 
radius  at  each  end,  the  wire  radius  being  assumed  to  vary  linearly 
along  the  wire.  The  gap  regions  are  assumed  to  be  conical  trans¬ 
mission  lines,  the  source  or  the  sink  being  located  at  the  apex  of  the 
cones.  The  wire  radius  must  be  defined  as  zero  where  a  gap  is  desired. 

Two  types  of  gaps  can  be  defined,  see  Fig,  1,  The  gap  at  one 
end  of  wire  section  1  connects  a  conical  wire  to  a  ground  plane,  while 
the  gaps  at  the  junction  of  wire  sections  5  and  6  connect  two  conical 
wires.  Any  number  of  wires  may  be  joined,  and  a  gap  be  defined  at 
each  wire  going  to  the  junction.  In  fact,  junctions  of  wires  are  handled 
most  accurately  by  the  program  when  a  short-circuited  gap  is  defined 
at  each  wire  of  the  junction.  In  each  gap  an  emf  is  defined  in  series 
with  a  generator  impedance  which  can  be  a  series-or  parallel-R-L-C 
circuit,  see  Fig.  1.  The  gaps  can  be  connected  by  any  number  of  four- 
terminal  networks,  see  Fig.  2.  The  program  could  handle  any  type  of 
linear  network  but  is  presently  set  up  to  handle  only  a  two-wire  trans¬ 
mission  line  (or  single  wire  if  the  network  is  ground-connected  as  in 
Fig.  2B). 

The  gaps  are  invariably  needed  in  antenna  problems.  Ir,  scatter¬ 
ing  problems,  the  source  is  not  an  emf  placed  in  the  wire  structure  but 
rather  a  distant  antenna;  e.g.,  a  simple  dipole  that  will  produce  a  plane 
wave  at  the  wire  structure.  In  order  to  enable  the  solution  of  scattering 
problems  as  well  as  antenna  problems,  dipole  current  sources  can  be 
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defined  anywhere  with  respect  to  the  wire  structure*  Each  dipole  source 
is  assumed  to  have  a  current  distribution  varying  as  Cos  (kx)  along  the 
dipole,  x  *  0  being  in  the  center  of  the  dipole,  and  k  being  the  free  space 
propagation  constant.  For  some  applications  it  may  even  be  desirable 
to  use  a  combination  of  gap-emf*s  and  dipole  current  sources. 

The  order  in  which  the  coordinates  of  the  ends  of  the  linear  wire 
sections  are  read  is  not  arbitrary.  The  reading  of  coordinates  must 
start  at  a  junction  or  an  open  end  of  the  wire  structure  and  must  pro¬ 
ceed  along  a  wire  until  the  next  junction  or  open  end  is  met.  Figure  1 
shows  one  possible  order  in  which  the  wire  sections  of  a  specific  antenna 
can  be  arranged.  The  positive  direction  of  current  is  assumed  to  be 
the  direction  in  which  the  wire  coordinates  are  read. 

In  many  applications  the  current  distribution  on  an  antenna  is 
known  to  possess  certain  symmetry  properties  equivalent  to  the  exist¬ 
ence  of  electric  or  magnetic  image  planes.  While  the  mutual  coupling 
effect  between  all  parts  of  the  antenna  must  be  taken  exactly  into  account, 
it  is  not  necessary  to  let  the  computer  solve  for  the  image  currents. 

Since  a  substantial  amount  of  computer  time  can  be  saved  by  making  use 
of  the  symmetry  properties  of  the  antenna,  thv  program  is  written  so 
that  the  antenna  can  be  imaged  in  the  planes  Z  *  0  and  Y  *0,  see  Fig, 3, 
The  plane  Z  *  0  can  be  defined  only  as  an  electric  "ground  plane,  " 
while  the  plane  Y  *  0  can  be  defined  as  an  electric  or  a  magnetic  plane. 

The  symmetrically  excised  circular  ring  array  is  another  example 

of  an  antenna  which  possesses  symmetry  properties  that  the  program 

will  utilize  at  a  considerable  saving  in  computer  time.  The  circular 

ring  array  is  assumed  to  have  N  antennas,  the  nrth  antenna  element 

being  constructed  by  rotation  of  the  first  antenna  element,  the  reference 

360 

element,  through  an  angle  —  (n-  1}  degrees  about  the  z-axis,  see 
Fig,  4,  A  ring  array  of  this  type  can  be  excited  in  N  linearly  indepen¬ 
dent  "azimuthal"  modes.  The  m*th  azimuthal  mode  is  characterized 
by  a  gap  voltage  in  the  n*th  antenna  element  that  is  related  to  the 
gap  voltage  in  the  reference  element  as  follows: 


Y  *  V,ej  N 
n  1 

With  this  excitation  of  the  ring  array,  the  current  distribution  In(s} 
in  the  n*th  element  will  be  related  to  the  current  distribution  I^(s)  in 
the  reference  element  as  follows: 

•  2ir  i  i  \ 

r  ,  .  -  ,  .  1  - rr—  in-  1  )m 

I  (s)  *  I,  (s)eJ  N  ' 

n  1 

where  s  is  a  coordinate  along  the  wire. 

The  symmetry  properties  of  a  symmetrically  excited  ring  array 
can  also  bi  made  use  of  to  great  advantage  when  the  excitation  of  the 
ring  array  is  arbitrary,  because  any  excitation  of  the  ring  array  can 
be  represented  by  a  linear  combination  of  the  symmetrical  modet  of 
the  ring  array. 

Three  types  of  azimuth  variations  of  the  ring  array  can  be 
handled  by  the  program.  The  exponential  variation  described  above 
is  the  most  general  variation  which  can  be  applied  to  any  ring  array 
according  to  Floquet*s  theorem.  The  program  will  also  handle  the 
cosine-  or  the  sine-variations  which  may  often  be  used  with  advantage 
when  such  variations  are  known  to  be  existing. 

i'he  program  "WIRTNA"  is  actually  written  for  a  ring  array. 

This  does  not  complicate  the  treatment  of  single  antennas  which  are 
merely  interpreted  as  ring  arrays  with  one  element. 

The  program  will  calculate  the  radiation  pattern  gain  in  ?ny  num¬ 
ber  of  directions.  The  direction  *  5  '--b^ervation  is  described  in  terms 
of  the  elevation  angle  9  and  the  azimuth  angle  <j>  as  shown  in  Fig.  4. 
Since  the  far-field  is  generally  elliptically  polarized,  the  gain  of  the 
major  axis  field  component,  the  gain  of  the  minor  axis  field  component, 
the  gain  of  the  vertical  field  component,  and  the  gain  of  the  horizontal 
field  component  are  computed  in  addition  to  the  total  gain.  The  axial 
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ratio  of  the  polarization  ellipse,  the  angle  between  the  major  axis  cf  the 
ellipse  and  the  horizon,  and  the  direction  of  rotation  of  the  field  vector 
are  listed  in  the  pattern  outputc  The  effect  of  a  finite  conductivity  of 
the  ground  can  be  taken  into  account  only  in  radiation  pattern  calcula¬ 
tions. 

The  program  provides  for  omitting  the  contributions  to  the  radia¬ 
tion  pattern  arising  from  the  dipole  current  sources.  The  omission  of 
these  contributions  to  the  radiation  pattern  is  usually  desired  in 
scattering  calculations. 

The  program  will  read  a  maximum  of  9  "pattern  cards,  "  each 
pa.ttern  card  defining  a  mesh  of  directions  of  observations  at  which 
the  pattern  gain  is  to  be  evaluated. 

The  program  will  handle  any  number  of  cases.  The  input  data 
cards  for  the  cases  to  be  analyzed  are  stacked,  and  the  last  case  is 
terminated  by  an  END-card. 
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The  input  data  for  the  prog  am  can  be  read  according  to  the 
following  sequence  of  FORTRAN  statements: 

DIMENSION  ITITLE(5) 

COMPLEX  RV<j)LT,  RAMP 

110  REAP  1,  (ITITLEUI),  1  m  1,  5)  IAZIM,  ISYMTRY,  NCURS<j)UR 
IF  (ITITLE«3HEND)  120t  200,  120 

120  READ  2,  NELEM,  NWIRE,  MSTART,  MMAX,  EPSGR,  C<|>NDGR, 
C<J)NDWR 

NWIRNCUR  *  Nw'IRE  +  NCURS(()UR 
D(}>  130  I  *  1, NWIRNCUR 

130  READ  3,  X1(I),  Y1(I),  Z1(I),  RADIUS  1(1),  X2j(I),  Y2(l),  Z2(I),  RADIUS  2(1} 
D(J)  190  NFR  *  1,  NFREQ 
READ  4,  FREQ 

READ  5,  (INTEVAL(I),  I  *  1,  NWIRE) 

D<t>  140  I  *  1,  NPAT 

140  READ  6,  IPATTYP(I),  NELEVA(I),  ELEVSTAR(I),  ELEVSTEP(I). 

NPHI(I),  PHISTAR(I),  PHISTEP(I),  ICURPATfl) 

D<t>  180  IEXC  x  1,NEXCITA 
D(j>  150  1 ,  NGAPEXC 

150  READ  7, 1,  RV$j/T(I),  YV<f>LT<I),  ISERPAR(I),  <t>HM(I),  HENRYM(I), 
FARADP(I) 

D<t>  160  I  *  1,NCURS#)UR 
160  READ  8,  RAMP(I) 

D<}>  170  I  *  1,  NETS 

170  READ  9,  ITYPENET(I)j(NETGAP(I,  J),  J  *  1, 4),  PARI  (I),  PAR2(I), 
PAR3  (I),  PAR4(1),  PAR5(I),  PAR6(I) 

180  CONTINUE 

READ  10,  LEND 
190  CONTINUE 
G<J>  T(t>  110 
200  CONTINUE 
END 


Format  Statements 


1  FORMAT  (5A8,  2  (AS,  2X)  15), 

2  F0RMAT  (415,  F10.0,  F10.4,  F10.0) 

3.  Fd>RiAA.T  (8F10.4) 

4.  FORMAT  (IX,  F9.4) 

5.  FORMAT  (1615) 

6.  FORMAT  (A8,  2X,  2(15,  2F10.4),  15} 

7.  FORMAT  (15,  C(F10.4,  F10.4)  A8,  2X,  3F10.4) 

8.  FORMAT  (IX,  C(F9.4,  F10.4)) 

9.  FORMAT  (A8,  2X,  412,  2X,  6F10.2) 

10.  FORMAT  (A3) 

While  the  above  sequence  of  FORTRAN  statements  is  a  valid 

method  of  reading  the  input  information,  it  should  also  be  pointed  out  that 
the  program  permits  a  substantial  saving  in  the  input  data  card  prepara- 
tion  when  calculations  of  an  antenna  are  to  be  made  at  several  frequencies 
for  the  same  number  of  integration  intervals,  the  same  patterns,  and  for 
the  same  excitation  modes.  In  such  cases  the  information  on  the  integra¬ 
tion  intervals,  the  patterns  and  the  excitation  modes  is  read  only  for  the 
first  frequency.  For  each  following  frequency  only  the  frequency  card 
need  be  read,  an^  the  program  will  automatically  assume  that  calcula¬ 
tions  are  to  be  made  for  the  same  number  of  integration  intervals,  the 
same  patterns,  and  the  same  excitation  modes  as  at  the  previous  fre¬ 
quency. 


An  explanation  of  the  symbols  used  in  the  above  sequence  of  FORTRAN 
statements  follows: 


C(f)NDGR 


is  the  conductivity  of  the  ground  z  *  0  over  which  the 
antenna  is  placed,  measured  in  mhos  per  meter. 

See  also  under  EPSGR. 
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C<(>NDWR 


is  the  conductivity  of  the  type  of  wire  composing  the 
antenna,  measured  in  mhos  per  meter.  If  C(|)NDWR 
is  left  blank,  the  wire  conductivity  is  assumed  to  be 
infinite. 

ELEVSTAR  is  the  first  elevation  angle  for  which  pattern  cal¬ 

culations  are  made. 

jst 

ELEVSTEP  is  the  difference  between  the  (n+1)  and  the  n*th 

elevation  angle  for  which  pattern  calculations  are 
made. 

EPSGR  is  the  relative  dielectric  constant  of  the  ground  z  *  0 

over  which  the  antenna  is  placed,  if  EPSGR  *  0,  the 
ground  is  assumed  to  be  perfectly  conducting.  If 
EPSGR  *  1,  the  antenna  is  assumed  to  be  placed  in 
free  space.  Otherwise,  the  ground  is  assumed  to  be 
perfectly  conducting  in  the  calculation  of  the  current 
distribution,  but  to  have  the  relative  dielectric  con¬ 
stant  EPSGR  in  the  radiation  pattern  calculations. 

FARADP(I)  is  the  capacitance  (C)  of  the  generator  impedance, 

measured  in  picofarads,  in  the  I*th  gap  of  the  refer¬ 
ence  antenna. 

FREQ  is  the  frequency  in  MHz.  The  card  from  which  FREQ 

is  read  must  have  an  "F"  in  the  first  column. 

HENRYM{I)  is  the  self  inductance  (L)  of  the  generator  impedance, 

measured  in  microhenrys,  in  the  I*th  gap  of  the 
reference  antenna. 
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IAZIM 


/ 


ICURPAT 


LEND 


is  the  azimuth  variation  assumed  in  the  ring  array- 
considered  and  may  be  any  of  the  following  symbols: 
EXPON  if  the  azimuth  variation  is  exponential 

COS  if  the  azimuth  variation  is  cosinoidal 

SIN  if  the  azimuth  variation  is  sinusoidal 

The  symbol  for  IAZIM  must  be  left  justified  in  the 

field  allowed.  If  a  single  antenna  is  considered,  the 
field  for  IAZIM  is  merely  left  blank. 


is  used  to  define  whether  the  contributions  from  the 


dipole  current  sources  to  the  radiation  pattern  are 


to  be  included  or  not,  ICURPAT  can  have  two  values: 


r 


ICURPAT  * 


if  the  contributions  from  the  dipole 
current  sources  are  to  be  included 
in  the  pattern. 

if  these  contributions  are  to  be  ex¬ 
cluded  from  the  pattern. 


*  3HEND. 


INTEVAL(I)  is  the  number  of  integration  intervals  used  on  the  I*th 

linear  wire  section  of  the  reference  element.  The  total 
number  of  intervals  between  one  junction  or  open  end 
and  another  junction  or  open  end  of  the  antenna  must  be 
equal  to  or  greater  than  3.  The  number  of  integration 
intervals  must  be  large  enough  to  assure  that  no  mte- 

,v  gration  interval  exceeds  about  0,2  wavelengths, 

* 

i 

\ 

t 

i 

i 
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IPATTYP 


ISERPAR(I) 


ISYMTRY 


'{ 


*  TITTLE 


is  used  to  define  various  forms  of  the  pattern,  gain 

output.  IPATTYP  can  be  any  of  the  following  symbols: 

DB  if  the  pattern  gain  in  dB  relative  to  an 

isotropic  antenna  is  desired. 

POWER  if  the  power  gain  relative  to  an  isotropic 
antenna  is  desired. 

VOLTAGE  if  the  amplitude  gain  relative  to  an  iso¬ 
tropic  antenna  is  desired. 

is  used  to  determine  whether  the  L,  C,  and  R  of  the 
generator  impedance  in  the  I*th  gap  oi  the  reference 
antenna  are  in  series  or  in  parallel  (see  Fig.  1). 

ISERPAR  can  be  any  of  the  follov  :ig  symbols : 

SERIES  if  a  series  combination  is  considered. 

PARALLEL  if  a  parallel  combination  is  considered. 

If  the  field  for  ISERPAR  is  left  blank,  the  generator 
impedance  is  assumed  to  be  0  independent  of  the 
values  read  for  L,  C,  and  R,  The  symbol  for 
ISERPAR  must  be  left  justified  in  the  field  allowed. 

is  used  to  define  the  presence  of  an  image  of  the  refer¬ 
ence  antenna  element  in  the  XZ-plane.  ISYMTRY  must 
be  one  of  the  following  symbols: 

ELECTRIC  if  the  imaging  assumes  the  XZ-plane  to  be 

an  electric  plane  {no  tangential  electric  field). 
MAGNETIC  if  the  imaging  assumes  the  XZ-plane  to  be  a 
magnetic  plane  (no  tangential  magnetic  field). 
If  the  field  allowed  for  ISYMTRY  is  left  blank,  no  image 
of  the  reference  element  is  introduced  in  the  XZ-plane. 

is  read  from  the  first  card  of  each  case  and  serves  as 
identification  for  the  customer.  The  information  on  this 
card  is  also  printed  on  the  output  listing. 
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«  ID 


rawssSEWBfl 


ITYPENET(I) 


MSTART 

MMAX 

NCURS(j)UR 

NELEM 

WELEVA 

NETGAP(I,  J) 


is  used  to  define  the  network  type  of  the  I*th  four-ter- 
minal  network  connecting  the  gaps  of  the  antenna.  In 
the  present  version  of  the  program,  ITYPENET  can 
only  be  the  following  symbol: 

TRLINE  which  stands  for  a  two-wire  transmission 
line  and  must  be  left  justified  in  the  field 
allowed. 


is  the  lowest  mode  number  of  the  ring  array  for  which 
calculations  are  to  be  made.  When  a  single  antenna  is 
considered,  choose  MSTART  =  0. 

is  the  highest  mode  number  of  the  ring  array  for  which 
calculations  are  to  be  made.  When  a  single  antenna  is 
considered,  choose  MMAX  a  0. 

is  the  number  of  dipole  current  sources  used. 

is  the  number  of  elements  considered  in  the  ring  array. 
When  a  single  antenna  is  considered,  choose  NELEM  *  1, 

is  the  number  of  elevation  angles  for  which  pattern  cal¬ 
culations  are  made.  The  gain  is  computed  NPHI  times 
for  each  value  of  the  elevation  angle. 

is  the  number  of  the  gap  to  which  the  J*th  terminal  of 
the  I?th  four-terminal  network  is  connected,  NETGAP 
can  have  values  from  1  through  NGAP,  but  can  also  have 
the  value  0,  which  is  used  to  indicate  a  ground  connection. 
The  ordering  of  the  four  terminals  of  a  network  is  as 
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NETGAP(I,  J) 
continued 


shown  in  Fig,  2.  For  the  example  in  Fig.  2a,  the 
four  values  of  NETGAP  would  be  4,  3,  6,  and  5.  For 
the  example  in  Fig.  2b,  the  four  values  of  NETGAP 
would  be  1,  0,  2,  and  0, 

A  180  degree  phase  shift  (physically  obtained  '  .y  a 
phase  inverter,  or  by  twisting  the  transmissic  n  line) 
can  be  introduced  by  interchanging  two  of  the  NETGAP 
values.  An  additional  180  degree  phase  shift  intro¬ 
duced  in  the  examples  in  Fig.  2a  and  Fig.  2b  could  be 
obtained  by  choosing  the  following  values  for  NETGAP: 

3,  4,  6,  and  5  for  example  2a;  and  1,  0,  0  and  2  for 
example  2b, 

NETS  is  the  number  of  four-terminal  networks  connected  to 

the  gaps  of  the  antenna.  This  number  is  not  an  input 
parameter.  The  program  keeps  automatically  a  count 
of  the  number  of  "network"  cards  read. 

NEXCITA  is  the  number  of  "excitation  modes"  for  which  calcula¬ 

tions  are  to  be  made.  Each  excitation  mode  corresponds 
to  a  particular  set  of  applied  emf’s  and  generator  im¬ 
pedances  in  the  gaps  of  the  reference  element  and,  if 
present,  a  particular  set  of  dipole  current  sources.  For 
a  single  antenna  with  only  one  gap,  NEXCITA  would  most 
often  be  1.  This  number  is  not  an  input  parameter.  The 
program  keeps  automatically  a  count  of  the  number  of 
excitation  modes, 

NFREQ  is  the  number  ox  frequencies  for  which  the  antenna  is  to 

be  analyzed.  This  number  is  not  an  input  parameter. 

The  program  keeps  automatically  a  count  of  the  number 
of  frequencies. 
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NGAP 


NGAPEXC 


NPAT 


NPHI 


NWIRE 


PARI  (I) 


is  the  number  of  gaps  on  the  reference  antenna  element. 
This  number  is  determined  by  the  computer,  and  the 
gaps  are  numbered  from  1  to  NGAP  in  the  order  in 
which  the  gaps  are  encountered  when  the  wire  co¬ 
ordinates  are  being  read.  The  double  gap  shown  in 
Fig,  1  counts  as  two  gaps, 

is  the  number  of  gaps  of  the  antenna  that  are  not 
short-circuited. 

is  the  number  of  "pattern"  cards  read.  This  number  is 
not  an  input  parameter.  The  program  keeps  automati¬ 
cally  a  count  of  the  pattern  cards  as  they  are  read. 

is  the  number  of  azimuth  angles  for  which  pattern  cal¬ 
culations  are  made.  The  gain  is  calculated  NELEVA 
times  for  each  value  of  the  azimuth  angle, 

is  the  number  of  linear  wire  sections  composing  the 
reference  antenna  element. 

is  the  resistance  (R)  of  the  generator  impedance, 
measured  in  ohms,  in  the  Pth  gap  of  the  reference 
antenna. 

is  a  quantity  that  describes  a  certain  parameter  of 
the  I*th  network  connecting  the  gaps.  For  the  trans- 
mission-line  network,  PARI  is  the  impedance  of  the 
transmission  line  (in  ohms). 
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PAR2(I) 


PHIS  TAR 


PHISTEP 


RADIUS  1(1) 


RADIUS  2(1) 


RAMP(I) 


RV<j)LT(I) 


XI  (I) 


is  a  quantity  that  describes  a  certain  parameter  of 
the  Pth  network  connecting  the  gaps.  For  the  trans¬ 
mission  line  network,  PAR2  is  the  length  of  the  trans¬ 
mission  line  (in  inches  if  all  other  dimensions  of  the 
antenna  were  read  in  inches  -  otherwise  this  length 
is  measured  in  wavelengths), 

is  the  first  azimuth  angle  for  which  pattern  calculations 
are  made, 

is  the  difference  between  the  (n+l)st  and  the  nTth  azi¬ 
muth  angle  for  which  pattern  calculations  are  made. 

is  the  wire  radius  at  the  first  end  point  of  the  Pth  linear 
wire  section  of  the  reference  antenna  element.  RADIUS 
1  is  measured  in  inches,  RADIUS  1  must  be  0  where  a 
gap  is  located. 

is  the  wire  radius  at  the  second  end  point  of  the  I*th 
linear  wire  section  of  the  reference  antenna  element. 
RADIUS  2  is  measured  in  inches.  RADIUS  2  must  be 
0  where  a  gap  is  located. 

is  the  complex  current  flowing  in  the  center  of  the  I’th 
dipole  current  source  and  is  measured  in  amperes. 

is  the  complex  emf  applied  in  the  ITth  gap  of  the  refer¬ 
ence,  antenna  and  is  measured  in  volts. 

is  the  x-coordinate  of  the  first  end  point  of  the  Pth 
linear  wire  section  of  the  reference  antenna  element. 

XI  is  measured  in  inches. 


X2(I)  is  the  wire  radius  at  the  second  end  point  of  the  I*th 

linear  wire  section  of  the  reference  antenna  element. 
X2  is  measured  in  inches. 

Y 1  (I)  is  the  y-coordinate  of  the  first  end  point  of  the  Ilth 

linear  wire  section  of  the  reference  antenna  element. 
Y1  is  measured  in  inches. 

Y2(I)  is  the  wire  radius  at  the  second  end  point  of  the  I*th 

linear  wire  section  of  the  reference  antenna  element. 
Y2  is  measured  in  inches. 

Z1(I)  is  the  z-coordinate  of  the  first  end  point  of  the  I*th 

linear  wire  section  of  the  reference  antenna  element. 
Z1  is  measured  in  inches. 

Z2(I)  is  the  wire  radius  at  the  second  end  point  of  the  I*th 

linear  wire  section  of  the  reference  antenna  element. 
Z2  is  measured  in  inches. 


Computer  Printout 


The  computer  program  will  calculate  the  current  distribution  on 
the  antenna,  as  well  as  the  antenna  input  impedance  and  the  radiation 
pattern.  The  computer  first  lists  the  input  data.  For  each  frequency  and 
each  excitation  mode,  the  computer  then  lists  the  gap  sources  and  the 
networks  connecting  the  gaps  (if  any).  For  each  azimuthal  mode  of  the 
ring  array,  the  computer  prints  the  amplitude  (in  amperes)  and  the 
phase  (in  degrees)  of  the  current  at  the  center  and  at  the  end  points  of 
each  integration  interval.  The  input  impedance  data  printed  consists  of 
the  input  impedance  and  admittance  seen  at  each  gap.  The  radiation 
pattern  printout  gives  information  on  the  absolute  gain  of  the  antenna  in 
the  directions  specified.  In  general,  the  radiation  field  will  be  elliptically 
polarized,  and  the  computer  prints  the  gain  of  the  vertical  field  com~ 
ponent,  the  gain  of  the  horizontal  field  component,  the  gain  of  the  field 
component  along  the  major  and  minor  axes  of  the  polarization  ellipse, 
as  well  as  the  total  gain.  The  computer  also  prints  the  axial  ratio  of  the 
polarization  ellipse,  the  angle  between  the  major  axis  and  the  horizontal 
plane,  and  the  sense  of  rotation  of  the  polarization  ellipse. 

The  computer  finally  prints  a  brief  summary  of  the  calculations. 
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INSTRUCTIONS  FOR  USE  OF  PROGRAM  "PLOTPAT" 

(A  Supporting  Program  for  Program  "WIRTNA") 

Control  Data  3800  Version 

Program  PLOTPAT  is  used  to  plot  antenna  patterns  using  the 
pattern  data  generated  by  the  linear  antenna  program  WIRTNA.  The 
pattern  data  are  punched  on  data  cards  in  a  very  compact  form  by 
WIRTNA.  These  cards  are  read  by  a  subroutine  BYNIN  in  program 
PLOTPAT  along  with  other  data  cards  which  indicate  the  particular 
pattern  to  be  plotted.  PLOTPAT  is  written  quite  generally,  and  any 
azimuth  or  elevation  pattern  cut  may  be  chosen.  Program  PLOTPAT 
executes  on  the  CDC  3800  and  places  the  pattern  plotting  instructions 
on  a  magnetic  tape.  The  instructions  on  this  magnetic  tape  are  then 
read  by  the  CDC  160  A  computer  which  drives  the  California  Computer 
Products  plotter  (CDC  Model  165). 


DESCRIPTION  OF  CARD  INPUT 


400  ID  s  0 

IW  =  0 

20  READ  I,  (TITLE(I),I  a  1,  10),  M(j)DE 
IF  (TITLE(l)-lH  )  30,100,30 

30  ID  =  ID+  1 

READ  2,  IGROUND,  EREF(ID),  NANG(1 ,  ID),  ANGSTAR(1,  ID), 
ANGSTEP{1,  ID),  NANG(2,  ID),  ANGSTAR(2,  ID),  ANGSTEP(2,  ID) 

CALL  BTNIN{  ,  P(IW+1 ),  NW<1)RDS(ID) 

IW  *  IW+1 
GO  TO  20 


C  ONE  BLANK  IS  READ  AFTER  BINARY  DATA  DECK  READ  BY  BYNIN 
IN  ORDER  TO  GET  TO  STATEMENT  100 

100  READ  3,  B<j)TT(j)M,  NCIRCLE 

110  READ  4,  TITL,  PATTYPE,  ANGLE,  BNDANGL,  BNDANGU, 
ITYPE,  IC<t>NTINU 

C  AN  UNLIMITED  NUMBER  OF  TITL  PATTYPE  CARDS  MAY  EE  READ  IN 
SERIES 

C  IN  ORDER  TO  GET  FROM  STATEMENT  110  TO  100,  READ  6HNEWB(J»T 
INTO  TITLE. 

C  IN  ORDER  TO  END  ALL  PLOTTING,  READ  3HEND  INTO  TITLE. 

C  IN  ORDER  TO  READ  NEW  PATTERN  DATA  FROM  ANOTHER  BINARY  DECK 
READ  A  BLANK  INTO  TITL. 


FORMAT  STATEMENTS 

1  FORMAT  (9A8,A4, 14) 

2  FORMAT  (  15,  El  5.  5,  2(15,  2F10.5)) 

3  FORMAT  (F10.  4,  3110) 

FORMAT  (A8,2X,A2,8X,3F10.4,2I10) 
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An.  explanation  of  the  symbols  used  in  these  statements  follows: 

TITLE  is  read  from  the  first  card  of  each  binary  data  deck  to 

identify  the  case. 

M(J)DE  is  the  order  of  the  ring  array  mode. 

IGR(|)UND  If  IGR(j>QND  =  1,  the  ground  is  present;  if  IGR(j>UND  »  0, 

no  ground. 

EREF  a  normalization  factor,  gain  =  E/EREF. 

NANG(1,I)  number  of  elevation  angles. 

NANG(2,I)  number  of  azimuth  angles. 

NW<j)RDS(I)  4.  0*NANG{1 , 1)*NANG(2, 1) . 

ANGSTAR(1,I)  the  first  elevation  angle  for  which  pattern  data  is  punched 
on  the  data  deck  . 

ANGSTEP(l.I)  the  increment  in  elevation  angle. 

ANGSTAR(2, 1)  the  first  azimuth  angle  for  which  pattern  data  is  punched 
on  the  data  deck  . 

ANGSTEP(2,I)  the  increment  in  azimuth  angle. 

P  an  array  in  which  the  pattern  data  is  stored.  The  total 

.number  of  words  available  in  this  array  is  12  102  words. 

B(j)TT(j)M  the  lowest  DB  value  on  the  DB  plot;  otherwise  BOTTOM 

is  meaningless. 

NCIRCLE  »  1  then  plot  the  unity  gain  cir  cle. 

NCIRCLE  s  0  then  normalize  pattern  to  unity  and  do  not  plot  the  circle. 

PATT\PE  is  a  Hollerith  constant  determining  whether  an  elevation 

or  azimuth  cut  is  to  be  plotted. 
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PATTYPE  a  EL  then  plot  pattern  versus  elevation  at  an  azimuth 

defined  by  ---angle. 

PATTYPE  a  AZ  then  plot  pattern  versus  azimuth  at  an  elevation 

defined  by  ---angle. 

PLOT  PATTERNS  BETWEEN  BANANGL  AND  BANANGU  (DEG. } 

ITYPE  determines  the  type  of  pattern  plot. 

ITYPE  a  1  power  plot. 

ITYPE  =  2  DB  plot 

ITYPE  a  3  AMP  plot. 

ITYPE  a  4  AMP  plot  of  major  and  minor  axes  in  place  of  THETA 

and  PHI  comp. 

makes  it  possible  to  plot  several  curves  on  one  graph. 


IC(j)NTINUE 
IC<|)NTINU  a  1 


then  add  to  the  previous  plot 
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,  Integral  equations  for  the  current  distribution  on  an  arbi-  ' 
'trary  wire  antenna  have  been  programmed  for  a  digital  computer,  j 
;The  resulting  computer  program  will  calculate  the  current  distri¬ 
bution,  the  input  impedance,  and  the  radiation  pattern  of  a  wire 
[antenna  of  arbitrary  geometry.  Considered  as  a  transmitting  an¬ 
tenna,  the  antenna  is  excited  by  electromotive  forces  applied  in 
any  number  of  gaps  along  the  wires.  The  program  developed  will 
.'permit  these  gaps  to  be  connected  by  an  arbitrary  non-radiating 
[network,  and  it  will  take  into  account  resistive  and  reactive 
•loading  along  the  wires.  The  program  will  also  treat  the  antenn^ 
as  a  receiving  antenna,  if  desired,  and  then  uses  a  distant  dipo] 
as  the  source  of  the  induced  current  distribution.  The  numerical 
integral  equation  method  used  is  essentially  an  exact  method.  The 
principal  errors  in  the  calculations  arise  in  the  numerical  inte¬ 
gration  of  the  integrals  in  the  integral  equations,  and  these 
errors  can  be  controlled.  It  is  usually  quite  easy  to  obtain 
accuracies  far  tetter  than  normal  measurement  accuracies,  and  at 
a  much  lower  cost. 
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13.  ABSTRACT  (cont) 


Numerical  results  are  presented  for  a  simple  dipole,  a  Yagi 
antenna  array,  a  log-periodic  antenna  array,  an  airplane  tail-cap 
antenna,  a  center-fed  whip  antenna,  and  a  compact  VHF  antenna. 
Most  results  have  been  verified  by  measurements. 
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